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ABSTRACT 


This  report  describes  the  results  of  a theoretical  study 
of  the  problem  of  a through  crack  In  a tensile  plate  speci- 
men. Experimental  observations  Indicate  that  the  crack  front 
grows  stably  in  the  specimen  Interior  forming  a "thumb  nail" 
shape  prior  to  unstable  fracture.  The  amount  of  interior 
growth  before  fracture  is  found  to  depend  on  both  geometry 
(relative  specimen  thickness),  the  material  properties,  and 
the  amount  of  yielding  which  tends  to  delay  the  onset  of  rapid 
fracture. 

The  goal  of  this  work  Is  to  study  the  interaction  of  ma- 
terial nonlinearity  with  geometry  and  their  combined  effect 
on  the  fracture  process.  To  this  end,  two-dimensional  elastic 
plastic  and  three-dimensional  elastic  calculations  have  been 
performed.  The  finite  element  method,  specialized  to  crack 
problems.  Is  the  tool  used  to  obtain  numerical  results.  Pre- 
dictions on  the  increments  of  growth  along  the  crack  front 
are  made  on  the  basis  of  the  strain  energy  density  theory 
that  assumes  crack  trajectory  to  coincide  with  path  of  mini- 
mum strain  energy  density  function.  Related  three-dimensional 
elastic-plastic  calculations  are  currently  being  developed  to 
complete  the  modeling  of  the  ductile  fracture  process  In 
plate  specimens. 
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TWO-DIMENSIONAL  ELASTIC-PLASTIC  ANALYSIS  OF  LONS  CRACKS 


A.  Introduction 

The  behavior  of  a center  crack  specimen  containing  a 
long  crack  and  correspondingly  short  ligament  Is  considered. 
Experimental  observations  of  crack  growth  through  the  last 
ligament  to  the  specimen  surface  Indicate  that  the  crack  turns 
away  from  Its  Initial  direction  and  final  failure  occurs  at 
an  angle  to  the  original  crack  plane.  Figure  1.  This  change 
of  direction  may  be  either  clockwise  or  counterclockwise  and 
the  choice  Is  believed  to  depend  on  a small  nonal Ignment  of 
the  crack  with  the  mid-plane  of  loading  not  treated  In  the 
analysis . 

A similar  phenomenon  Is  observed  In  the  fracture  of  mod- 
erately thick  specimens  where  the  combinations  of  material 
properties  lyleld  stress)  and  specimen  geometry  violate  the 
ASTM  requirements  for  Klc  testing  [1].  Here,  the  crack  grows 
In  the  plane  of  the  Initial  crack  configuration  in  the  speci- 
men Interior  but  deviates  from  that  plane  near  the  specimen 
surfaces  forming  "shear  lips".  A three-dimensional  descrip- 
tion of  the  fracture  sequence  Is  pertinent.  The  crack  grows 
stably  from  Its  initial  configuration  In  a planar  fashion 
developing  a curved  crack  front  as  shown  In  Figure  2.  The 
point  of  Instability  and  the  breakthrough  of  the  curved  crack 
front  to  the  plate  surfaces  (at  an  angle  to  the  plane  of  the 
main  crack)  appear  to  occur  at  about  the  same  point  In  the 

loading  history.  Again,  local  conditions  are  assumed  to  de- 
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terralne  whether  the  shear  Up  will  form  on  one  crack  surface 
or  the  other. 

The  authors  hypothesize  that  crack  turning  upon  approach- 
ing a specimen  surface  In  the  cases  of  (a)  the  "last"  liga- 
ment of  fracture  and  (b)  the  formation  of  shear  lips  are  dif- 
ferent manifestations  of  the  same  phenomenon.  Further,  the 
understanding  of  the  process  Involved  In  the  formation  of 
shear  lips  Is  fundamental  to  the  prediction  of  the  stable- 
unstable  transition  In  fracture.  That  problem  Is  of  consid- 
erable Interest  to  the  authors;  however.  It  is  difficult  and 
complex  to  manage  because  It  requires  three-dimensional  elas- 
tic-plastic analysis  for  a specimen  containing  a growing  crack. 
Therefore,  the  phenomenon  of  crack  turning  during  the  "last" 
ligament  of  growth  Is  studied  here  using  two-dimensional  anal- 
ysis to  gain  Insight  Into  the  shear  lip  formation  process. 

B.  Stress  Analysis 

A two-dimensional  finite  element  computer  program  developed 
In  [2,3]  Is  employed  to  carry  out  the  stress  analysis.  The 
program  Is  based  on  the  twelve  node  Isoparametric  element. 

It  Includes  a specialized  crack  tip  element  which  enforces 
the  Hutchlnson-RIce  plastic  singularity  solution  [4].  The 
singularity  solution,  and  therefore  the  finite  element  pro- 
gram, Is  based  on  a power  hardening  model  of  material  behavior 
and  the  J2  deformation  theory  of  plasticity,  l.e.,  the  strain- 
stress  relationship  Is  expressed  In  the  form 
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where  e^j  3 e^j  - ^ epp5ij  the  stra*n  deviation  and  the 

stress  deviation  is  s^j  = ^ appl51j‘  The  Hutchlnson- 

Rlce  plane  strain  solution  for  the  asymptotic  stress,  strain, 
and  displacement  (relative  to  the  crack  tip)  fields  respec- 
tively takes  the  form: 


'ij  * °yd  Kr  "?T  S1J<e> 


E1J  * r °yd 


n n n 1 


n n " 1 

ui  " I ayd  K r Q1(8) 


where  (r,8)  represent  a polar  coordinate  system  centered  at 
the  crack  tip  and  where  the  plastic  intensity  factor  K de- 
pends on  global  geometry  and  loading  conditions. 

The  crack  tip  element  Is  semi-circular  In  shape  and  Its 
displacement  shape  function  Is 
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where  u and  K are  the  generalized  nodal  displacement  compo- 
o „ 
nent  and  K£  * Kn. 

C.  Results 

The  center  crack  specimen  Is  considered  with  two  width  to 
length  ratios:  4:2  as  employed  In  the  previous  report  [5]  and 

4:4.  Table  1 shows  the  crack  lengths  and  applied  loads,  oQ, 
which  were  treated. 


TABLE  1 - APPLIED  STRESS  TO  YIELD  STRESS  RATIO  FOR 
DIFFERENT  SIZE  SPECIMENS 
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The  values  of  the  material  parameters  chosen  for  these 
calculations  are: 

E = 1.0,  v = 0.3 

CTyd  = a * 0*5  and  n * 5.0 

Figure  3 shows  a uniaxial  tensile  stress-strain  curve  for  this 
material  model.  Typical  grid  patterns  for  the  two  aspect  ra- 
tios are  shown  In  Figures  4 and  5.  Figures  6 to  22  show  the 
regions  in  each  specimen  which  have  yielded,  i.e.,  the  circled 
dots  represent  quadrature  points  at  which  the  applied  stress 
Is  greater  than  the  yield  stress.  Also  Included  In  these  fig- 
ures Is  a curved  line  from  the  crack  tip  to  the  specimen  bound- 
ary. This  line  represents  the  numerical  result  for  the  minimum 
strain  energy  density  path  along  which  the  crack  Is  assumed  to 
fol low. 

D.  Discussion 

In  carrying  out  these  computations,  the  authors  were  par- 
ticularly interested  in  obtaining  information  about  the  plastic 
region  between  the  crack  tip  and  the  specimen  boundary  In  order 
to  get  Insight  Into  the  shear  lip  problem.  The  results  obtained 
Indicate  that  In  every  case  there  exists  an  elastic  zone  between 
the  crack  tip  and  the  boundary.  However,  for  the  width  to 
length  ratio  of  4:2,  a secondary  plastic  zone  is  observed  to 
form  at  the  free  boundary.  This  secondary  plastic  zone  Is  ab- 
sent from  the  results  with  specimen  width  to  length  ratio  of 

-6- 


4:4.  These  results  emphasize  the  Importance  of  geometry  on 
stress  field  and  therefore  crack  propagation. 

The  pertinent  question  Is  whether  the  crack  would  run 
through  the  region  of  maximum  plastic  deformation  toward  the 
specimen  boundary  or  If  It  would  attempt  to  avoid  the  deformed 
material,  growing  Into  the  elastic  zone  In  front  of  the  tip  and 
then  diverging  from  the  line  of  symmetry  to  break  through  the 
last  ligament  at  an  angle.  The  paths  of  minimum  strain  energy 
density  consistently  lie  between  the  region  of  maximum  plastic 
deformation  and  the  horizontal  axis  (connecting  the  crack  tip 
to  the  free  surface).  Thus,  the  strain  energy  density  theory 
predicts  that  the  crack  will  propagate  close  to  the  Inner  elas- 
tic-plastic boundary  Intersecting  the  surface  at  an  oblique 
angle.  Careful  experimental  measurements  are  needed  to  con- 
firm the  accuracy  of  this  prediction. 

Three-dimensional  elastic-plastic  analysis  for  a specimen 
with  a curved  crack  front  Is  underway  to  examine  the  similarity 
of  the  shear  lip  breakthrough  phenomenon  with  the  two-dimen- 
sional "last  ligament"  problem. 


THREE-DIMENSIONAL  ELASTIC  ANALYSIS 


AND  GROWTH  PREDICTIONS 

A three-dimensional  element  which  Incorporates  the  elastic 
crack  front  singularity  and  Is  applicable  to  curved,  as  well 
as  straight,  crack  fronts  Is  formulated.  A number  of  calcula- 
tions are  performed  to  establish  the  accuracy  of  results  which 
can  be  anticipated  using  the  specialized  elements  In  conjunc- 
tion with  isoparametric  elements.  Then  the  finite  element 
procedure  Is  employed  together  with  a strain  energy  density 
based  criterion  to  predict  Increments  of  crack  growth. 

A.  Formulation  of  a Specialized  Finite  Element  Procedure  for 

Three-Dimensional  Crack  Problems 

The  finite  element  displacement  method  Incorporating  Iso- 
parametric elements  [6]  will  be  employed.  In  general,  an  ele- 
ment Is  enclosed  by  six  surfaces  with  nodal  points  chosen  along 
the  element  edges  (a  five  surface  element  will  be  Introduced 
later  In  the  text).  A coordinate  transformation  Is  employed 
to  map  the  element  Into  a cube.  Figure  23.  The  Inverse  mapping 
function  can  be  expressed  In  the  form 


■ m 
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y 

* £ Ni (p,s ,t) 
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where  1 ■ 1,  the  number  of  nodes  associated  with  the  element, 

and  (x^.y^.z^)  are  the  coordinates  of  the  1th  node.  The  func- 
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tlons  N^Cp.s.t)  associated  with  the  mapping  are  polynomial 
functions  with  the  property 


1 for  p,s,t 


Nj  (p»s,t) 


[o  for  p,s,t 


p^.Sy.t*  (the  point  in  the 

mapped  shape  which 
corresponds  to  x. , 

yi * z1> 

(5) 


j * 1 


The  same  functions  N^(p,s,t)  are  used  In  the  approxima- 
tion for  the  displacement  field  within  the  element,  l.e.,  the 
Interior  displacement  components  are  expressed  In  terms  of 
the  corresponding  nodal  components  as 


“ • 
ux 

uy 

* £ Nj(p,s»t) 

X1 

■ uz. 

L“z,J 
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These  functions  N^(p,s,t)  are  termed  as  "displacement  shape 
functions" . 

The  strains  components  are  obtained  In  terms  of  the  nodal 
displacement  components  by  differentiating  equation  (6),  i.e.. 
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where  the  Jacobian  of  the  transformation 
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is  inverted  to  obtain  the  needed  partial  derivatives.  The  re- 
sult is  expressed  symbolically  as 


m t 
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where 


e - (ex,  ey » ez.  Yxy,  Yyz,  Yzx) 
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The  stress  components  are  related  to  the  strain  components 
through  the  generalized  Hooke's  law 
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The  strain  energy  in  an  element  can  then  be  expressed  In  the 
form 


W « \ $T  [k]  t (12) 

where  [k]  Is  the  element  stiffness  matrix  given  by 

[B]T  [D][B]  detJ  dpdsdt  (13) 

and  where  detJ  is  the  determinate  of  the  Jacobian  for  the 
transformations  given  in  equation  (8). 

The  potential  energy  of  the  body,  calculated  as  the  sum  of 
the  element  strain  energy  contributions  minus  the  work  done  by 
applied  forces.  Is  minimized  with  respect  to  nodal  displace- 
ment components  to  obtain  a set  of  linear  algebraic  equations 
for  the  determination  of  their  values. 

Specialized  "singular"  elements  are  developed  to  properly 
model  the  crack  front  singularity.  The  approximate  form  for 
the  displacement  field  within  elements  adjacent  to  the  crack 
edge  Is  chosen  to  Include  the  first  term  of  the  asymptotic 
solution 
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ux  = kj(z)/2fr  [ (5-8v ) cos  | - cos  + 0(r) 

= k1 (z)  ux(r,e)  + 0 (r) 

uy  * ^4^  k-j  (z)/2r  [ ( 7 - 8 v ) sin  | - sin  + 0 (r ) (14) 

= k-j  (z)  Uy(r,0)  + 0 (r) 

uz  = 0(r)  as  r + 0 

where  E Is  Young's  modulus  and  v Is  Poisson's  ratio. 

Here,  (r,8,z)  are  local  cylindrical  coordinates  centered 
on  the  crack  front  with  z tangent  to  the  crack  edge.  The 
geometry  of  these  singular  elements  is  designed  to  follow  the 
local  cylindrical  coordinate  system  in  which  the  singular  so- 
lution is  given.  A typical  singular  element  configuration  Is 
shown  in  Figure  24.  Its  cross  section  in  the  plane  normal  to 
the  crack  edge  is  pie  shaped.  These  singular  elements  are 
then  arranged  to  fill  a volume  with  circular  cross  section 
centered  on  the  crack  edge.  Figure  25. 

The  assumed  form  for  the  displacement  field  within  each  of 
the  singular  elements  consists  of  two  parts.  It  contains  the 
standard  Isoparametric  element  approximation  to  the  displace- 
ment field  plus  a contribution  from  the  elastic  singular  so- 
lution. In  order  that  this  displacement  field  matches  that 
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of  standard  Isoparametric  elements  at  nodes  common  to  both 
elements,  the  contribution  from  the  singularity  solution  Is 
forced  to  vanish  at  these  nodes.  The  resulting  displacement 
field  can  be  written  using  local  coordinates  (x,y,z)  along  the 
crack  front  In  the  form 

ux(p,$,t)  3 l Ni(p,s,t)ux^  + k1 (z)[ux(r,6) 

- I uy(r1 ,01 )N1 (p.s.t)] 

uyCp»s,t)  = l N.j  lp»s ,t)u  + k]  (z)[u  (r,e) 

1 (15) 

“ I 3y(ri,ei)N1(p,s,t)] 

uz(p,s,t)  3 l N1(p,s,t)uz 

where  uz(r,0)  and  Oy(r,0)  are  given  by  equation  (14). 

The  variation  of  the  stress  Intensity  factor  along  the 
crack  front  Is  approximated  In  a piece-wise  manner  consistent 
with  the  finite  element  displacement  approximations,  l.e., 

•<!  (z ) * ^ N*(s)ki  (16) 

where  s Is  the  variable  In  the  mapped  space  which  corresponds 
to  z and  where  the  summation  Is  over  nodes  adjacent  to  the 

ft 

element  along  the  crack  edge.  The  shape  functions  N^ts)  are 
given  by  N*(s)  « N^-l.s.-l). 
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Transforming  to  the  global  coordinate  system  (X,Y,Z),  the 
assumed  form  for  the  displacement  field  in  a singular  element 
becomes 


Nj(p,s,t)ux  + | {N*(s)[ux(r,e)cosa 
Nj(p,s,t)ux(ri ,6i )cosai]}ki 
Njfp.s.Oiiy  + | (N*(s)[Uy(r,0) 
N1(p,s,t)uy(rj,6j)])ki 
NitP.s,t)uz^  + J (N*(s)[uxCr,0)s1na 
Nj  Cp»s,t)ux(rj ,6j  JsInajlJk^ 


07) 


where  a Is  the  angle  between  the  normal  to  the  crack  edge  (In 
the  plane  of  the  crack  front)  and  the  x axis  given  by 
a = arc  tan  (-  fff|f) . 

The  strain  components  are  obtained  by  differentiating  the 
displacement  expressions  In  equation  (17)  as: 
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This  approximation  to  the  strain  field,  used  In  the  elements 
adjacent  to  the  crack  edge,  contains  the  elastic  strain 
singularity  plus  non-singular  contributions  associated  with 
standard  isoparametric  elements.  The  nodal  values  of  the 
stress  Intensity  factors  can  now  be  treated  as  generalized 
nodal  displacement  components  or  degrees  of  freedom  to  be 
determined  through  the  finite  element  procedure.  Let  $ be 
a vector  representation  of  the  degrees  of  freedom  associated 
with  a singular  element,  then 

^ ~ , Uy  , u ^ , Uy  ,.«.,  u_  , k, , k9» . . . ,k_ 

A1  rl  Z1  x2  zm  Z n 

where  m Is  the  number  of  nodes  of  the  element  and  n is  the 
number  of  these  nodes  which  lie  along  the  crack  edge.  Then 
the  approximation  to  the  strain  components  in  a singular  ele- 
ment can  be  expressed  In  the  same  form  as  that  employed  for 
a standard  element,  i.e., 

e - £8*]  t*  (19) 


and  the  corresponding  element  stiffness  matrix  Is  given  as 


[k] 


[D]  £B*J  detJ  dpdsdt 


(20) 


B.  Remarks  on  Accuracy  - Straight  Crack  Front  Results 

The  specified  finite  element  procedure  described  above  has 
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been  employed  to  solve  the  problem  of  a center  crack  (straight 
crack  front)  specimen  under  uniform  tensile  loading.  Results 
were  obtained  for  two  specimen  thicknesses  In  order  to  study 
the  Influence  of  specimen  thickness  on  the  through-the-thlck- 
ness  variation  of  the  near  edge  fields.  Initial  results  were 
based  on  the  grid  patterns  shown  In  Figure  25.  (The  crack 
edge  was  taken  as  straight  with  a singular  element  radius  of 

0. 25  or  1/1 6th  of  the  half  crack  length.  The  corresponding 

curves  for  the  variation  of  the  normalized  stress  Intensity 
k -i 

factor  ( ) with  the  normalized  transverse  coordinate  (z/h) 

a /a 

are  g1ven°1n  Figure  26.  The  results  shown  In  this  figure  are 
for  the  case  where  the  stress  Intensity  factor  Is  constrained 
to  be  zero  at  the  free  surface.  Note  the  expected  trend, 

1. e.,  as  the  plate  thickness  h is  decreased  relative  to  the 
in-plane  dimensions,  the  degree  of  variation  of  the  stress 
Intensity  factor  Increases.  (For  an  Infinitely  thick  plate, 
the  stress  Intensity  factor  Is  constant  at  the  plane  strain 
value). 


These  calculations  have  been  repeated  using  a signifi- 
cantly smaller  singular  element  radius,  l.e.,  rn  = 0.05  or 

kl 

1/40  of  the  half  crack  length.  The  curves  for  against 

a / a 

z/h  are  shown  In  Figure  27.  Note  that  less  variation  of 
stress  intensity  with  distance  from  the  center  plane  of  the 
plate  Is  predicted  for  both  plate  thicknesses,  and  further, 
that  the  results  do  not  change  significantly  from  one  plate 
thickness  to  the  next.  The  average  Interior  value  for  the 
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stress  Intensity  factor  Is  also  Influenced  by  the  change  in 
grid  pattern. 

Figure  28  shows  the  Influence  of  the  surface  constraint 
on  the  variation  of  the  stress  Intensity  factor  across  the 
plate  thickness.  Note  that  this  Influence  Is  only  local, 
l.e.,  near  the  plate  surface.  The  extent  of  the  observed 
difference  between  the  results  for  the  stress  Intensity  factor 
constrained  and  free  at  the  plate  surface  Is  Influenced  by 
grid  pattern.  As  the  thickness  of  grid  layers  near  the  sur- 
face Is  decreased,  the  constraint  is  found  to  have  a more 
local  effect. 

These  results  give  the  reader  a measure  of  the  Inaccura- 
cies In  predictions  which  can  result  even  with  the  use  of 
singularity  elements.  The  explanation  for  these  discrepancies 
Is  related  to  the  use  of  the  singularity  solution  to  model  the 
solution  throughout  the  singular  elements.  The  asymptotic  so- 
lution, equation  (14),  Is  applicable  in  regions  close  to  the 
crack  edge,  l.e.,  distances  r which  are  small  compared  to  the 
other  geometric  dimensions.  For  the  first  grid  pattern,  the 
core  element  radius  Is  not  sufficiently  small  In  comparison 
with  the  plate  thickness,  and,  as  a result,  the  plane  strain 
condition  assumed  for  the  singular  solution  Is  violated  over 
a significant  portion  of  these  elements. 

These  results  demonstrate  that  the  region  over  which  the 
plane  strain  singularity  holds  Is  directly  related  to  plate 
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thickness.  Thus,  If  a crack  growth  criterion  is  applied  at 
some  absolute  distance  rQ  from  the  crack  edge,  increased 
variation  of  growth  along  the  crack  front  will  be  predicted 
as  plate  thickness  Is  decreased.  This  prediction  Is  consis- 
tent with  experimental  observations. 

C.  Crack  Growth  Modeling 

Calculations  have  been  performed  for  various  crack  front 
shapes  in  conjunction  with  crack  growth  modeling.  The  results 
Indicate  that  the  variation  of  the  stress  Intensity  factor 
along  the  crack  edge  Is  sensitive  to  the  details  of  the  crack 
front  geometry.  This  suggests  that  there  may  be  a "natural" 
crack  front  shape,  which  will  necessarily  depend  on  material 
properties,  specimen  geometry,  and  loading,  at  which  unstable 
fracture  can  occur.  The  concept  is  supported  by  experimental 
observations  [7]  for  specimens  of  only  moderate  thickness. 

When  a fatigue  crack  is  grown  from  an  Initial  straight  cut, 
the  crack  front  quickly  develops  a curved  shape  which  Is  then 
maintained  during  further  fatigue  crack  propagation.  The 
shape  of  the  crack  front  Is  influenced  by  the  amplitude  of 
the  loading  cycles  with  generally  increased  curvature  accom- 
panying Increased  load  amplitude.  If  a monotonic  loading  Is 
then  applied,  crack  growth  Initiates  stably  In  the  specimen 
Interior  thereby  Increasing  the  front  curvature.  Instability 
and  the  development  of  shear  lips  (nonplanar  crack  front 
growth)  are  observed  to  occur  at  about  the  same  time. 
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Consider  the  crack  growth  process  under  roonotonic  load- 
ing. The  load-load  point  displacement  record  for  a crack 
growth  process  Involving  a moderate  thickness  specimen  can  be 
characterized  by  a linear  portion,  a nonlinear  regime  and  an 
Instability  (Figure  29).  In  the  linear  range,  no  significant 
growth  or  yielding  is  occurring.  On  the  other  hand,  during 
the  nonlinear  (yet  stable)  portion  of  the  curve,  both  plastic 
deformation  and  stable  growth  may  be  occurring  simultaneously. 
The  relative  contributions  of  these  two  phenomenon  will  de- 
pend on  material  properties,  geometry  (thickness),  and  load. 

The  crack  growth  components  clearly  dominates  at  instability. 

It  appears  necessary  to  be  able  to  characterize  the  sepa- 
rate contributions  of  plastic  deformation  and  stable  crack 
growth  and  to  understand  how  they  vary  with  geometry  and  ma- 
terial properties.  In  order  to  be  able  to  predict  the  crack 
growth  process  in  a structural  component  based  on  specimen 
test  data. 

The  goals  of  this  research  effort  are  to  model  the  stable 
growth  process  described  above  by  separately  analyzing  the 
Influences  of  crack  front  geometry  (local  growth)  and  the 
contributions  from  material  non! Inearl ty . The  subsequent  aim 
is  to  account  for  these  two  contributing  factors  slmul taneously 
and  thereby  predict  the  crack  growth  process  under  monoton ic 
loading  In  terms  of  material  properties  and  geometry. 
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In  this  section,  results  of  three-dimensional  elastic 
finite  element  analyses  are  used  In  conjunction  with  a local 
growth  criterion  to  predict  increments  of  crack  growth  and 
correspondl ng  changes  in  crack  front  shape.  (Current  efforts 
Including  plastic  Influences  will  be  presented  In  a later 
report).  The  approach  employed  here  yields  meaningful  Infor- 
mation concerning  crack  front  shape  and  changes  in  specimen 
compliance;  however.  It  is  incapable  of  treating  the  question 
of  stability  because  It  does  not  account  for  material  irre- 
verslbl 1 1 ty . 

D.  Local  Crack  Growth  Criterion 

In  general,  both  the  direction  and  magnitude  of  local 
crack  growth  from  the  current  crack  front  are  expected  to 
vary  with  position  along  that  front.  The  strain  energy  den- 
sity field  surrounding  the  crack  front  is  the  basis  for  the 
development  of  the  local  crack  growth  criterion  to  be  con- 
sidered here  [8-10].  It  is  postulated  that  the  path  of  growth 
from  each  point  along  the  crack  edge  will  follow  the  minimum 
strain  energy  density  path  emanating  from  that  point.  Further, 
growth  at  points  along  the  current  crack  front  will  Initiate 

when  the  strain  energy  density  at  a "core"  distance  (rQ)  along 

* 

the  minimum  path  reaches  a preset  value  W . 

Hence,  the  continuous  growth  of  the  crack  front  Is  approx- 
imated by  discrete  Increments  of  growth.  The  amount  of  growth 
at  a point  along  the  crack  front  in  an  Increment  is  taken  to 
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Note,  a current  limitation  of  the  finite  element  crack 
tip  element  formulation  discussed  earlier  Is  that  the  crack 
front  must  be  normal  to  the  plate  surface  at  the  point  of 
Intersection.  (Other  three-dimensional  crack  tip  elements 
such  as  the  quarter  point  elements  [11]  and  Tracey's  pie 
shaped  elements  [12]  are  also  subject  to  this  limitation). 
Minor  adjustments  of  the  predicted  shape  of  crack  growth  In- 
crements are  imposed  to  satisfy  this  restriction. 

Some  additlon&l  Information  of  potential  Interest  obtained 
from  the  finite  element  results  includes: 

(a)  The  variations  of  the  stress  intensity  factor  along 
the  crack  front  at  each  increment  - see  Figure  35.  The  rapid 
variations  predicted  In  the  vicinity  of  the  free  surface  are 
believed  to  be  a result  of  forcing  a 90°  intersection  between 
the  crack  front  and  the  specimen  surface. 

(b)  The  variations  of  the  dlstortional  strain  energy  den- 
sity component  directly  ahead  of  the  existing  crack  front 

Cat  6 = 0°)  shown  In  Figure  36. 

The  fact  that  the  dlstortional  component  of  strain  energy  den- 
sity, normally  associated  with  plastic  deformation.  Increases 

★ 

as  the  free  surface  Is  approached  and  does  so  more  rapidly 


The  decrease  In  dlstortional  strain  energy  density,  In  the 
Immediate  vicinity  of  the  surface  may  be  the  result  of  numeri- 
cal Inaccuracies. 


i. 


mi 


for  curved  crack  fronts  than  straight  ones  Implies  the  develop- 
ment of  plastic  constraint  at  the  surface  and  the  formation  of 
shear  lips. 
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Figure  8 - Yield  pattern  for  2:1  aspe 
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Figure  9 - Yield  pattern  for  2:1  aspect  ratio  and 
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Figure  11  - Yield  pattern  for  1:1  aspe 
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Figure  15  - Yield  pattern  for  1:1  aspect  ratio  and 
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Figure  19  - Yield  pattern  for  1:1  aspect  ratio  and  a * 0.28  a ..  £•  = 0.9 


0 ‘02 


20.0 


Surface  of  plate 


Figure  26  - Variation  of  stress  intensity  factor  along  crack  front. 


Surface  of  plate 


The  Influence  of  the  surface  constraint  on  the 
stress  intensity  factor  variation. 
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Figure  32  - Strain  energy  density  variation  along  a curved  crack  front  - 
first  increment  at  crack  growth. 
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Figure  35  - Variation  of  stress  Intensity  factor 
along  crack  front  as  It  grows. 
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